Ginsburg-Landau theory of supersolid 
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We develop a simple Ginsburg-Landau theory to study all the possible phases and phase transi- 
tions in He, analyze the condition for the existence of the supersolid (SS) and map out its global 
phase diagram from a unified framework. If the condition favors the existence of the SS, we use the 
GL theory to address several experimental facts and also make some predictions that are amenable 
to experimental tests. A key prediction is that the X-ray scattering intensity from the SS ought to 
have an additional modulation over that of the NS. The modulation amplitude is proportional to 
the Non-Classical Rotational-Inertial (NCRI) observed in the torsional oscillator experiments. 



1. Introduction: A solid can not flow. While a 
superfluid can flow without any resistance. A supersolid 
(SS) is a new state of matter which has both the solid and 
superfluid order. The possibility of a supersolid phase in 
*iJe due to very large zero point quantum fluctuations 
was theoretically speculated in 1970 [1-4]. Over the last 
35 years, a number of experiments have been designed 
to search for the supersolid state without success. How- 
ever, recently, by using torsional oscillator measurement, 
a PSU group lead by Chan observed a marked 1 ~ 2% 
Non-Classical Rotational Inertial ( NCRI ) of solid *_ff e 
at ~ Q.2K, both when embedded in Vycor glass and in 
bulk "^He [5] . The NCRI is a low temperature reduction 
in the rotational moment of inertia due to the super- 
fluid component of solid ^He [3]. The authors suggested 
that the supersolid state is responsible for the NCRI. 
The PSU experiments rekindled extensive theoretical in- 
terests [6-8] in the still controversial supersolid phase of 
^He. There are two kinds of complementary theoretical 
approaches. The flrst is the microscopic numerical simu- 
lation [6]. The second is the phenomenological approach 
[7,8]. In this paper, by constructing a Ginsburg Landau ( 
GL ) theory, we will address the following two questions 
: (1) What is the condition for the existence of the SS 
state ? (2) If the SS exists, what are the properties of the 
supersolid to be tested by possible new experiments. We 
develop a simple GL theory to map out the ^i7e phase 
diagram and study all the phases and phase transitions in 
a unifled framework. If the repulsive coupling in Eqn.3 in 
the GL theory is sufficiently small, the SS phase becomes 
stable at low enough temperature. The resulting solid at 
high pressure is an incommensurate solid with zero point 
vacancies whose condensation leads to the formation of 
the SS. The theory can be used to address several phe- 
nomena observed in the PSU experiments and also make 
sharp predictions to be tested by possible future experi- 
ments, especially X-ray scattering experiments in the SS 
state. 

2. Ginsburg-Landau theory of ^ife: Let's start by 
reviewing all the known phases in ^iJe. The density of a 
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tice vector. In a normal liquid (NL), if the static liquid 
structure factor S{k) has its first maximum peak at fc„, 
then near fc„, S~^{k) ~ r„ -|- c(fc^ — k'^Y . If the liquid- 
solid transition is weakly first order, it is known that the 
classical free energy to describe the NL-NS transition is 
[10]: 
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where r^ 



Tn + c(G'^ — k^^Y is the tuning parameter 
controlled by the pressure. 

It was known that the Superfluid (SF) to Normal Liq- 
uid transition at finite temperature is a 3d XY transition 
described by: 



u = K\^i^\''+m''+u\i;\^ + 



(2) 



where ip is the complex order parameter and t is the tun- 
ing parameter controlled by the temperature. 

The coupling between n{x) and '!/'(^) consistent with 
all the symmetry can be written down as: 
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where n*s, 



where the interaction must be repulsive wi > 0. Due to 
the lack of particle- hole symmetry in the NS, additional 
terms like n{x)ip^dr'4' should exist and is very impor- 
tant at zero temperature and will be investigated in [11]. 
However, in the classical phase transitions investigated 
in this paper, this term can still be neglected. As shown 
in section 5, the repulsive interaction is proportional to 
the temperature shift Tsf — Tss in Fig.l. 

In an effective GL theory, n(x) and ip{x) emerge as 
two independent order parameters. The total density of 
the system is nt{x) — n{x) + \'ip{x)\'^ where n(x) is the 
normal density and |?/'(a;)|^ is the superfluid density. A 
NS is deflned hy n^ ^ 0,< ip >= 0, while a SS is deflned 
by "-(5 7^ 0, < "0 >¥" 0- From the NL side, one can 
approach both the NS and the SF. Inside the NL, i > 0, 
■0 has a gap, so can be integrated out, we recover the 
_g and G is any non-zero reciprocal lat- NS-NL transition tuned by r^ in Eqn.l ( Fig.l ). Inside 



normal solid (NS) is deflned as n(x) = Uq + J^g ^g^ 



the NL < n{x) >= no, so we can set ng = for G 7^ 
in Eqn.3, then we recover the NL to SF transition tuned 
by t in Eqn.2 ( Fig.l ). 

Although the NL-NS and NL-SF transitions are well 
understood, so far, the SF-NS transition has not been 
investigated seriously. This transition may be in a com- 
pletely different universality class than the NL-NS tran- 
sition, because both sides break two completely different 
symmetry: internal global (7(1) symmetry and transla- 
tional symmetry. It is possible that the solid reached 
from the SF side is a new kind of solid than the NS 
reached from the NL side. In the following, incorporat- 
ing quantum fluctuations into f^{ see Eqn.4 ) and con- 
sidering the repulsive coupling between ip and n sector 
in Eqn.3, we will determine the global phase diagram of 

3. Two-component Quantum GL theory in the 
Tp sector: In this section, we develop a two-component 
Quantum GL theory in the ip sector to replace Eqn.2 to 
describe the superfluid side oi^He. The superfluid is de- 
scribed by a complex order parameter ip whose conden- 
sation leads to the Landau's quasi-particles. Although 
the bare "^He atoms are strongly interacting, the Lan- 
dau's quasi-particles are weakly interacting. The disper- 
sion curve of superfluid state is shown in Fig.l b which 
has both a phonon sector and a roton sector. In order 
to focus on the low energy modes, we divide the spec- 
trum into two regimes: the low momenta regime fc < A 
where there are phonon excitations with linear dispersion 
and high momentum regime \k — kr\ < A <C k^ where 
there is a roton minimum at the roton surface k = kr. 
We separate the complex order parameters iI){x,t) = 
ipi{x,T) + ^p2{x,T) into Vi(^,t) = J^ -f0je^''-^i:{k,T) 

and ■02(x,t) = I\k^k,\<A-(i:;^e'''''^^p{k,T) which stand 
for low energy modes near the origin and kr respectively. 
For the notation simplicity, in the following, /^ means 

I\k-k |<A- '^^'^ QGL action in the ip sector in the {k,uj) 
space becomes: 
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In the q -^ Q hmit, S{q) ~ q, u}{q) ^ q recovers the ipi 
sector near q = 0. The first maximum peak in S{q) cor- 
responds to the roton minimum in to{q) in the 'ip2 sector, 
namely, fc„ — kr- As one increases the pressure p, the 
interaction u also gets bigger and bigger, the first maxi- 
mum peak of S{q) increases, the roton minimum A^ gets 
smaller and smaller. Across the critical pressure p = Pc, 
there are two possibilities (1) The resulting solid is a com- 
mensurate solid, then < t/j >= (2) The resulting solid 
is an in-commensurate solid with vacancies even at T = 
whose condensation leads to < ^ >^ [1,2,7]. Which 
case will happen depends on the strength of the repul- 
sive coupling vi and will be analyzed in the next section. 
Case (1) is trivial, the SS phase in Fig.l does not ex- 
ist. In the following, we only focus on case (2). From 
Eqn.l and Eqn.4, we can see that n and tp2 have very 
similar propagators, so the lattice formation in n sector 
with n{x) — ^g n^e^^''^ and the density wave forma- 
tion in ■ip2 sector with < -02(2?) >= e^^^ J2m=i AmC*'^™''^ 
where Qm = kr happen simultaneously. The "02 sector 
alone is described by n = 2 component {d + l,d) with 
d = 3 Lifshitz action [10]. The repulsive coupling in 
Eqn.3 f 1 > simply shifts the DW by suitable constants 
along the three unit vectors in the direct lattice. These 
constants will be determined in the next section for dif- 
ferent n lattices. Namely, the SS state consists of two 
inter-penetrating lattices formed by the n lattice and the 
ip superfluid density wave (SDW). 

Combining the roton condensation picture in this sec- 
tion with the results in section 2, we can sketch the fol- 
lowing global phase diagram of '^He in case (2). 
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where t ^ T ~ Tsf where Tsf ~ 2.17 K is the critical 
temperature of SF to NL transition at p = 0.05 bar and 
^r ^ Pc ^ P where p,, ^ 25 bar is the critical pressure of 
SF to the SS transition at T = 0. 

4. SF to SS transition and global phase dia- 
gram: In the SF state, the Feymann relation between 
the Landau quasi-particle dispersion relation in the ip sec- 
tor and the static structure factor in the n sector holds: 



Fig.l: (a) The phase diagram in case (2). T controls thermal 
fluctuations, while P tunes quantum fluctuations. The SS exists 
only when the repulsive coupling vi in Eqn.3 is sufliiciently small. 
Thick ( thin ) lines are 1st ( 2nd ) order transitions. The critical 
temperatures of NL to SF and NS to SS transitions drop slightly as 
the pressure p increases because of the quantum fluctuations [11]. 
(b) The separation of low ( phonon ) and high ( roton ) momenta 
regime in the SF. 



5. The NS to SS transition: In this section, we 
approach the SS phase from the NS side and determine 
its lattice structure. In the NL, the BEC happens in the 
■01 sector at /c = 0, -02 has a large gap and can be sim- 
ply integrated out. In the NS, 4) stands for the vacancies 
due to large zero point quantum fluctuations [1,2,7]. Due 
to the n lattice formation, the mass of ^"1 was increased 
hq = T — Tss- Subtracting it from 
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t + v^^no = T-TsF leads to TsF -Tss ^ rioAvi. 
For simplicity, one can set w^^ = 0. If vi is larger than 
a critical value v > Vci, then Tss is suppressed to zero. 
This is the case (1) discussed in the last section. The 
SS disappears from the phase diagram Fig.l. The re- 
sulting solid is a commensurate solid. In the following, 
we only focus on vi < Vd, so the SS state exists in the 
Fig.l with Tss > which is the case (2) in the last 
section. The resulting solid is an in-commensurate solid 
with vacancies even at T = whose condensation leads 
to the formation of the SS. Then the temperature shift 
AT = TsF — Tss is an effective experimental measure 
of the repulsive coupling vi in Eqn.3. In the presence of 
the periodic potential of n{x) lattice, V' will form a Bloch 
wave, the u self-interaction in the ip sector in Eqn.4 will 
favor extended Bloch wave over strongly localized Wan- 
nier state. In principle, a full energy band calculation 
incorporating the interaction u is necessary to get the en- 
ergy bands of ip. Fortunately, qualitatively physical pic- 
ture can be achieved without such a detailed energy band 
calculation. In the following, substituting the ansatz 
< Mx) >= ae'^^ and < V2(x) >= e*^^ E™=i A^e*'^™-? 
where Qm = Q into Eqn.3, we study the effects of n lat- 
tice on Ip ~ ipi + il;2- In order to get the lowest energy 
ground state, we must consider the following 4 condi- 
tions: (1) because any complex ?/; ( up to a global phase ) 
will lead to local supercurrents which is costy, we can also 
take Ip to be real, so Q„i have to be paired as anti-nodal 
points. P has to be even (2) as shown from the Feymann 
relation Eqn.5, Qm, m — 1, • • • , P are simply P shortest 
reciprocal lattice vectors, then translational symmetry 
of the lattice dictates that e{K = 0) = e{K ~ Qm), 
■01 and "02 have to condense at the same time. (3) The 
point group symmetry of the lattice dictates Am — A 
and is real (4) The repulsive interaction wi > favors 
-0(a; = 0) = 0, so the Superfluid Density Wave (SDW) 
p ~ 1-01^ can avoid the n lattice as much as possible. It 
turns out that the 4 conditions can fix the relative phase 
and magnitude of "01 and 'ip2 to be ^2 ~ Oi+tt, A = a/ P, 
namely, -0 = ae*^(l - p Z]m=i cosQm ■ x). In this state, 
the crystal momentum fc = and the Fourier compo- 
nents are ^^{K = 0) = a, il){K — Qm) = —a/P which 
oscillate in sign and decay in magnitude. In principle, 
higher Fourier components may also exist, but they de- 
cay very rapidly, so can be neglected without affecting 
the physics qualitatively. In the following, we will dis- 
cuss 4 common lattices sc, fee, bcc, hep respectively. 



(a) P — Q: Qi,i — 1,2,3,4,5,6 are the 6 shortest 
reciprocal lattice vectors generating a cubic lattice. The 
maxima of the DW ipmax = 2a appear exactly in the mid- 
dle of lattice points at the 8 points a = ^(±ai ±a2±a3). 
They form the dual lattice of the cubic lattice which is 
also a cubic lattice. 

(b) P ~ 8: Qi,i ~ l,---,8 form the 8 shortest re- 
ciprocal lattice vectors generating a bee reciprocal lattice 
which corresponds to a fee direct lattice. The field is 
ip{x) — a[l~j{coaQi-x+cosQ2-x+cosQs-x+cosQ4-x)], 
The local superfiuid density p^^ — \tp{x)\'^. The max- 
ima of the DW Tpmax — 2a appear in all the edge centers 
such as (1/2, 0, 0) etc. and the centers of any cube such 
as (1/2,1/2,1/2). 

(c) P = 12: Qi,i = I,---, 12 form the 12 shortest 
reciprocal lattice vectors generating a fee reciprocal lat- 
tice which corresponds to a bee direct lattice. The field 
is ip(x) = o[l — i(cosQi • X + cos (32 • a? -I- cosQs • x + 
cos (54 • X + cos Q^ ■ X + cos Qe ■ x)]. The maxima of the 
DW ipmax = 4/3a appear along any square surrounding 
the center of the cube such as (1/2, /3, 0) or (1/2, 0, 7) etc. 
Note that ^He in Vycor glass takes a bee lattice 

(d) Unfortunately, a spherical kr ~ Q surface can not 
lead to lattices with different lengths of primitive recip- 
rocal lattice vectors such as a hep lattice. This is sim- 
ilar to the classical liquid-solid transition described by 
Eqn.l where a single maximum peak in the static struc- 
ture factor can not lead to a hep lattice [10]. Here we 
can simply take the experimental fact that n forms a hep 
lattice without knowing how to produce such a lattice 
from a GL theory Eqn.l. Because for an idea hep lattice 
e/a — ■y/8/3, an hep lattice has 12 nearest neighbours, so 
its local environment may resemble that of an fee lattice. 
We expect the physics ( except the weak anisotropy of 
NCRI in the hep lattice discussed in [11] ) is qualitatively 
the same as that in fee direct lattice. 

Let's look at the prediction of our theory on X-ray scat- 
tering from the SS. For simplicity, we take the se lattice 
to explain the points. The other lattices will be discussed 
in [11]. For a lattice with j = 1, • • • , n basis located at dj, 
the geometrical structure factor at the reciprocal lattice 
vector K is S{K) = Ej=i ./j(^)e'-^'^"'^ where fj is the 
atomic structure factor of the basis at dj . The X-ray scat- 
tering amphtude I{K) ^ \S{K)\'^. For the SS in the se 
lattice, K = ^{Tiii + n2J + n^)k,di =0,^2 = ^(i+j + k), 
then S{K) = 1 + /(-i)"i+"2+n3 .^yj^ere / ~ p^"^ ~ o^. 
It is 1 -I- / for even K and 1 — / for odd K. In fact, due to 
large zero-point motion, any X-ray scattering amplitude 
I(K) will be diminished by a Debye-Waller factor. As 
shown in [11], the lattice phonon modes m in -0 are loeked 
to those of n, so there is a common Debye-Waller factor 
^ e^3^ <" > for both even and odd K. We conclude 
that the elastie X-ray scattering intensity from the SS 
has an additional modulation over that of the NS. The 



modulation amplitude is proportional to the maxima of 
the superfluid density p™"^ ^ a} which is the same as 
the NCRI observed in the PSU's torsional oscillator ex- 
periments. Unfortunately, so far, the X-ray scattering 
data is limited to high temperature T > 0.8if > T55 
[12]. X-ray scattering experiments on lower temperature 
T < Tss are needed to test this prediction. 

The results achieved in this section indeed confirm 
Fig.l achieved from the roton condensation picture in 
the last section. A low energy effective action involving 
the superfluid phonon 9, the lattice phonons u and their 
couplings will appear in a future publication. 

6. The vortices in the SS: In the sectors of ^, there 
are topological defects in the phase winding of 9 which 
are vortices. At T <C Tss, the vortices can only appear 
in tightly bound pairs. However, as T ^ ^551 the vor- 
tices start to become liberated, this process renders the 
total NCRI to vanish above T > Tss in the NS state. In 
the SF phase, a single vortex energy costs a lot of energy 
E^^ = ^ — ^InrM- where m is the mass of '^He atom, 
R is the system size and ^5^ ^ a is the core size of the 
vortex. This energy determines the critical velocity in SF 
v^^ > 30cm/s. In the SS state, because in the center of 
the SS vortex core, -0 = 0, so the vortices prefer to sit on 
a lattice site of the n lattice. Because the long distance 
behavior of SS is more or less the same as SF, we can es- 
timate its energy E^^ = ^ — ^In-^. We expect the core 

size of a supersolid vortex S,ss ^ 1/^ ^ 1/^r ~ a ~ £^sf- 
so inside the SS vortex core, we should also see the lattice 
structure of n. In fact, we expect that ^ss is of the order 
of the average spacing between the vacancies in the SS. It 
is interesting to see if neutron or light scattering experi- 
ments can test this prediction. Compared to E^^ , there 
are two reductions, one is the superfluid density, another 
is the increase of the vortex core size ^ss ^ £,SF- These 
two factors contribute to the very low critical velocity 
wf ■^ ~ 30/iTO/s. 

7. Discussions on PSU's experiments: Although 
the NS to SS transition is in the same universality class 
as the NL to SF one [8], it may have quite different off- 
critical behaviours due to the SDW structure in the SS 
state. We can estimates the critical regime of the NS to 
SS transition from the Ginsburg Criterion t^^ ^ ^|^AC 
where AC is the specific jump in the mean field theory. 
Because of the cubic dependence on £^ss, large £,ss leads 
to extremely narrow critical regime, the 3D XY criti- 
cal behavior is essentially irrelevant, instead mean field 
Gaussian theory should apply. This fact can be used to 
address two experimentally observable effects of the He3 
impurities. The first one was already observed by the 
PSU experiments and the second one is underway [13]. 
(1) The unbinding transition temperature Tss is deter- 
mined by the pinning due to impurities instead of by 
the logarithmic interactions between the vortices. So the 
'^He impurities effectively pin the vortices and raise the 



unbinding critical temperature Tss- On the other hand, 
'^He impurities will certainly decrease the superfiuid den- 
sity in both the ipi and ip2 sector just like^iJe impurities 
decreases superfluid density in the ^iJe superfluid. (2) 
The mean field specific heat jump at T = Tss may be 
smeared by the presence of ^He impurities. 

8. Conclusions: The GL theory developed in this pa- 
per leads to a global and unified picture of '^He physics 
at any temperature and pressure. If the repulsive cou- 
pling in Eqn.3 is sufficiently small, the SS state becomes 
a ground state at zero temperature. The solid at high 
pressure is an incommensurate solid with zero point fluc- 
tuations generated vacancies whose condensation leads 
to the formation of the SS. Assuming this is indeed the 
case, we investigate the SS state from both the SF and 
the NS side and find completely consistent description of 
the properties of the SS state. We found that the SF to 
the SS transition is a first order transition driven by the 
collapsing of roton minimum in the SF side, while the 
NS to SS transition is described by a 3d XY model with 
much narrower critical regime than the NL to SF tran- 
sition. The SS state is a uniform two-component phase 
consisting of a normal lattice plus a commensurate su- 
perfiuid density wave (SDW). The SDW in the SS state 
leads to a modulation on the X-ray scattering intensity 
over that of the NS. The modulation amplitude is pro- 
portional to the NCRI observed in the PSU's torsional 
oscillator experiments. 
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